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DISLOCATION AND BATE-PROCESS lEEOBIES 
By A. S. Nowlok and E. S. Maohlln 


SQMMABY 

(An equation for tlie steady-state rate of creep has heon derived, 
hy applying the theory of dislocations to the creep of pu3^e metals^ 
The form of this equation is In good agreement vlth empirical equa- 
tlons descrlhlng creep rates. The theory was also used to predict 
the dependence of steady-state rate of creep on physical constants 
of the material and good agreement vas obtained vlth data In the 
• literature for pure annealed metals, ^e rate of creep vas found to 

decrease vlth Increasing modulus of rlgldltj^ &hls fact sviggeets 
that one of the roqulrementa for a heat-resisting alloy Is that Its 
matrix be a metal that has a hl^ modulus of rigidity and therefore 
a hl^ modulus of elasticity.]^ 

DTERODDCTION 

The development of the gas turbine In combination vlth the jet 
as a power plant for military aircraft has focused attention on the 
need for heat-resisting alloys. ^One of the criteria used to evaluate 
heat-resisting alloys is creep reeletanoe’j) that Is^ the reelstaxioe 
to plastic deformation over a period of time. Cuiront evaluation of 
creep resistance Is accomplished by creep tests in which creep curves 
(elongation plotted egalnst time) are obtained at constant stress and 
temperature. A typical creep curve is shown In figure 1. The initial 
stage, in which the slope of the curve or rate of creep Is rapidly 
decreasing, 1s cooiaonly designated the primary or transient stage; 
secondary or steady-state creep refers to the stralght-llne portion 
of the creep cuarve. After a sufficient length of time, the rate of 
creep increases In the region designated the tertiary stage. In 
^ certain ceisos, the rate of creep continuously Increases emd cannot 

be divided into these three stages (reference 1(a) ) . The steady- 
state rate of creep is generally the criterion by which the creep 
^ resistance of heat-resisting alloys la expressed. Because this rate 

rapidly Increases with temperature, creep boocmes an Important factor 
In limiting safe-operating temperatures. 
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One of the first attempts to analyze steady-state creep by other 
than emplz*loal methods vas by Eanter (reference Z) in 1938. Further 
vork of a fundamental, nature vas carried out by E^zmann (refer- 
ence 3) vho applied Syrlng's theory of liquid flow, vhlch is a special 
application of the theory of rate processes, to the steady-state flow 
of metals. The problem of creep has also been attacked by the use of 
the theory of dislocations. A surrey of results obtained thus far by 
this approach is repozised by Seitz and Bead (reference 4) . Although 
fthese three theoretical treatments predict equations for the steady- 
state rate of creep that showed the correct dependence on stress end 
temperature, quant itat Ira predictions cannot be made from them as to 
how physical and structural constants of materials affect the creep 
rate*^ Those treatments therefore are not of groat value in leading 
to the synthosls of new materials that will have better creep prop- 
erties than those cuzrently used for heat-resisting purposes. 

In order to determine the dependence of steady-state creep rate 
on physical constants of materials, the investigation reported herein, 
which continuos a fundamental approach to the problem of creep, vas 
conducted at the NACA Cleveland laboratory in the early part of 1945. 
Chi equation for the steady-state rate of creep as a function of 
applied stroBS and temperature is derived by applying Byring* s theory 
of rate prooessos and the theory of dislocations to the problem^ 

This investigation is part of a program being conducted at the Cleve- 
land laboratory to evaluate the physical prppertlos of heat-reslstlhg 
alloys in tozms of physical constants, which are either known or 
easily measurable, in ozder to minimize the number of tests as well 
as to make possible the prediction of new compoeitlons and structures 
for better heat-resisting alloys than those currently used. 


SIMBOL8 

The following symbols are used in the theoretical analysis; 
d]_ distance between atoms in slip direction 

dg Interplanar spacing (fig. 4(a)) 

B modulus of elasticity 

F(T) function of temperature and material 

AFq^ free energy of activation per molecule 

f fraction ■vdiose value is about l/Z 

G modulus of rigidity at any temperature 
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H 




Gq modulus of rlgldltT* at absolute zero 
h Planck's constant, 6.62 x 10"®^ erg seconds 
heat of aotlTatlon per molecule 

k Boltzmann's constant, 1.38 x 10~ ° ergs per molecule per 

L distance betveen Is^exf actions In a single crystal (of the 
order of 1 micron) 

N^ number of dislocations per unit area that intersect a plane 
noxmal to the slip plane and containing the slip direction 

Ng number of sources of dislocations per unit volume 

p probability of occurrence of oscillation In orystallogra^lc 
dix^ctlon under consideration 

?■*" ratio of sources of dislocations at vhlch a positive activated 
complex has formed to total number of right-hand sources of 
dislocations; that is, probability of formatlcm of a positive 
activated complex at rl^t-hand source 

P" probability of formation of a negative activated complex at 
right-hand source 

P» - -k logg p^ 

q, stress concentration factor 

r rate of reaction per unit concentration of reactants 

B rate of generation (number of dislocations generated/seo for 

^ single generating source) 

Bg rate of generation of positive dislocations at rl^t-hond 
source 

B~ rate cf generation of negative dislocations at rlght-heind 
® source 

entropy cf activation per molecule 
t avereige time for activated complex to pass Into product form 

T absolute temperature, ^ 
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u tensile creep rate at steadj-state conditions 
u' shear rate 

T average velocity cf motion of dislocations 

T volume associated vlth one atom 

W elastic energy 

energy to be attained by thexmal motion 
X ratio of d^ to dg 

a temperature coefficient of modulus of rigidity 

& displacement shovn In figure 5(a) 

H Poisson's ratio 

a externally applied tensile stress 

T ^ear stress 

T'b back stress 

Te externally applied shear stress 

Tg average shear stress at sources of dislocations 

shear stress to be attained by thermal motions 
Tq critical local shear stress for generation of dislocation 

THEOET 

The features of the theory of rate processes and the theory of 
dislocations required for the development of the creep equation are 
first reviewed; the derivation of the equation for the steady -state 
rate of creep then follows. 

Eyrlng's Theory of Bate Processes 

The theory of rate processes (inference- 5} , os developed by 
Eyrlng and others, considers a reaction or any rate process to be 
the result of the crossing of a potential -energy barrier by molecules 
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whose energies have exceeded a certain minlnoim value. The reaction 
consists of the formation of an “activated coai^>lex“ capable of cross- 
ing the barrier followed by the passage of this complex over the 
barrier. (See fig. 2 .) The most important assumption Is that the 
initial reactants and activated complexes are always in egulllbrlum. 
Application of thermodynamic and statistical mechanical considera- 
tions shows that the number of activated complexes passing over the 
barrier per second (or rate of reaction) per imlt concentration of 
reactants is given by 

r . (1) 

The term kr/h can be re^urfed as the effective frequency at 
which activated complexes cross over the barrier; the exponential 
factor represents the probability of formation of an activated 
complex. In the calculation of this probability term, l^o contri- 
bution due to the translational degree of freedom along the "reaction 
coordinate" (the most favorable reaction path on the potential-energy 
surface) has been disregarded because it is included in the factor 
kr/h. The free energy of activation is Interpreted as an 

ordinary free -energy term and can be expressed by 

AF^ - AHa - TASq (2) 

where the heat of activation per molecule AE^^ la the hel^t of the 
potential-energy barrier in figure 2. 


General Features of Dislocation Theory In Bolatlon to Slip 

The theory of dislocations arose in relation to the mechanism 
of slip in single crystals of metal. Slip is the plastic defoxmatlcxi 
that occurs in an arbitrary short period of time when a metal is 
subjected to stress. This deformation is clutraoterlzed by the dis- 
placement of one part of the crystal lattice relative to another 
along a particular crystallographic plane and in a definite crystal- 
lographic direction. These slip pianos and slip directions are 
detected by the occurrence of slip bands, step-like discontinuities 
on the surface of single crystals that have been subjected to stress. 
The spaclngs of the bends are generally of the order of 1 micron 
(reference 1(b)). For any given crystallographic plane, there is a 
more or less critical value of the component of the applied stress 
above which the rate of plastic deformetioa rapidly increases. This 
value is called the critical shearing stress. Sarly attempts to 
explain the mechanism of slip based upon the assumption that the 
process consists in the sliding of pcu’allel planes over one another 
lead to theoretically determined values of the critical shearing 
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B tresses, uhloh vere cf the order of lOGO times too large. It there 
fore hecame eTldent that the true nature of slip Involves deviations 
from a perfect lattice. 

Of several mechanisms suggested the cm& that has been most 
successful In explaining the slip process Is the theory cf dislo- 
cations. This thecoy, •which was advanced by G-. I. Taylor (refer- 
ence 6), Orowan (reference 7), and Polanyl (reference 8), proposes 
that local deviations from a perfect lattice called dlslooatlone 
exist In single crystals, and -that movement of dlslocatltxis -bhrough 
the stressed orys-tal produces slip. A dlslooatlcHi consists In a 
stable arrangement of atoms such -that. In a region of a few atomic 
distances, n -i- 1 atvxos in the slip direction face n atoms across 
the slip plane. The point at •which the atoms are one-half on atomic 
spacing "out of step" Is called -the center of -the dislocation. The 
extension of -the dislocation In the slip plane and normal to the 
slip direction Is called the length of the dislocation. 

A simplified model of a dislocation In i^loh two neighboring 
planes are shown In cross section Is presented In figure 3; the 
plane of the figure Is normal to the length of the dislocation. A 
small shear stress applied as shown in figure 3(a) will cause a 
displacement cf the center of -the dislocation by one atomic distance 
(fig. 3(b)). The dislocation finally passes ccn 5 )letely out cf -the 
crystal and a perfect lattice Is restored, differing from the orig- 
inal lattice in that one plane lias now been translated relative to 
•the other by one atomic dls^tmco. In the process Just described, 
the entire linear dislocation must move as a unit. 

Barrett (reference 9) presents a sketch token from an earlier 
work by Taylor (reference 6) showing the -fcwo types of dislocation 
that can exist (fig. 4). If, In the nel^borhood cf -the center of 
the dislocation, more atoms are above the slip plane, a positive 
dislocation occurs and. If more atoms are below, a negative dislo- 
cation occurs. Fot the same shear stress, posltl^ve and negative 
dislocations move In opposite directions but the net result of •bhe 
passage -bhrou^ -the crystal Is 'the same, a translation of one part 
of the lattice -with respect to "the o^fchor by one atomic distance 
(fig. 4). 

Inasmuch as 'the peieslng of a dislocation out of -the crystal 
lattice results In •the loss of the dlslooatlcsi, now dislocations 
must be generated in order for continuous plastic deformation to be 
possible. A dynamic theory of dlelocations is therefore noceseary 
and •was de-ve loped by Orowan (reference 7) and extended by others. 
Including Seitz and Bead (reference 4). Because of the great energy 
req,ulred, the process of generation generally takes place at regions 
of hl^ stress oonoentratlon resulting f rem Imperfections In the 
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single orjstal.s. The ezlstenoe of such Imperfections as a result 
of accidents of growth during the fonnablon of single czystc^s Is 
generalljr accepted, hut the nature of the flaws that exist Is still 
a matter of oontrorersy. tar example, the theory that single crys- 
tals are mad.e up of a mosaic block structure Is d.iscussed. In refer- 
ence 4. Z-ray evidence Indicates that the average spacing of the 
Imperfection Is of the ord.er a£ 1 micron. Tl:e general dynamic 
picture obtained. Is ae follows: Dislocations are generated under 

shear stress at the regions of hi^ stress concentration, which will 
be called sources cf dislocations. A source of dislocations will be 
called a right-hand source If It generates dislocations to the ri^it 
and a left-hand source If it generates dislocations to iiie left. 

The dislocations eis they arc generated ere only a few atoms long 
because of the extreme low probability of generating a full-length 
dislocation (one whose length Is of the order of the spacing butween 
Imperfections) (reference 4) . As the dislocations move throu^ the 
crystal under shear sti-ess, they Increase In length and move in the 
slip plane and slip direction until they become stuck at an Imper- 
fection or In the vicinity of other dislocations. 

Fuarbher properties of dislocations and their use to explain the 
qualitative phenomena of plastic dof ormatlcai ere summarized In a 
aeries of four papers by Seitz and Bead (reference 4) , 


Derivation of an Equation for Creep Bate 

In the analy'sls of the creep process, it will be assumed as In 
reference 4 that creep os well as slip takes place by the motion of 
dislocations. The dynamic model previously described will bo used. 

The procedure la as follows: An equation will be derived for the 

steeidy-state rate of creep, which will show that the rate of genera- 
tion of dislocations la the prlmaiy factor that determines the varia- 
tion of creep properties among matorlt^s. ihe rate of generation 
will then be calculated In terms of rate-process theory as a frequency 
factor multiplied by the probability that a source of generation will 
attain the energy required to generate a dislocation. 

For a given block of volume Ip, If It is assumed that no 
’’annihilation" (union of pairs -jf dislocations of opposite sign) 
takes place and that steady-state creep is being dealt with, the 
number of dislocations generated por second within the block must 
equal the number leaving the block per second. The number of dis- 
locations leaving the block por second Is the number of dislocations 
In an area vL so 
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Inasmuch as the displacement produced vhen a dislocation passes out 
of the block is d^, the shear rate u* is given by 

u' = i (N^vDdj^ - N^vdi (4) 

The occurrence cf slip bands at distances apart of the order of 
1 micron Indicates an average cf one source cf dislocations for each 
imperfection; that is 

Ng - lA.^ (5) 

When equations (3), (4), and (S) are oooiblned, the expression for 
shear rate becomes 

U- - ^ Bg (6a) 


Up to this point, only the creep of single crystals has been 
considered. In a coisideration of polycrystalline metals, it vould 
appear that two types of creep may take place; creep occurring 
within the Individual grains or single crystals, and creep resulting 
from, intergranular motion. As noted in reference 4, present data on 
the effect of grain size on creep leaves doubt as to the possibility 
at the occurrence of intergranular creep. For example, the work 
reported in reference 10 on copper cf various grain sizes shows no 
detectable difference in creep strength other than that resulting 
from oxidation. Unfortunately, this work was carried out only at a 
single temperature. Creep in polyciystalllne metals is therefore 
assumed to be predominately the result of defcamnation occurring 
within the individual grains by the motiem. cf dislocations. If the 
orientation of the grains is such that slip occurs in the piano of 
maximum shear stress, then, as in iTOferenoo 3, the strain rata u' 
can be corrected to the tensile creep rate u by an additional 
factor cf about two-thirds. Such a factor would bo of no consequence 
here because d^/L is known only in ardor of magiitude; accurate 
values cannot be given until more is known about the nature of the 
crystal imperfections. Equation (te) can therefore be written 


approximately as 



(6b) 


The rate of generation E_ will now be calculated, A specific 
model for the mechanism of generation of a dislocation must be set 
up and certain simplifying asstm^tiens made. The approximations that 
will be made are similar to those used by Orowan (reference 7) but 
the treatment presented heire is based upc»i a more definite model them 
he used. Fu3rthermore, considerations in refezence 7 were not applied 
to steady-state creep. 
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The generated dislocation will be taten to te 1 ataa long. 

(The length of a dislocation Is Its extension in the slip plane and 
perpendicular to the slip direction.) The process of generation 
■will involve displacements of atoms lying along the slip direction 
by large fractions of an atomic distance. A simple model for the 
generation of a poslti-ve dislocation at a rl^t-hand source is pre- 
sented in flgu3?ee 5(a) and 5(b). If atoms 1 and 2 are s-tmul-baneously 
displaced In the directions shown in figure 5(a) by large fractions 
of the lattice spacing and atoms 3, 4^ and 5 are moving in the slip 
direction as f^own, the generation of a dislocation becomes possible. 
Before a dislocation ■will form, however, a roarrangement of the 
nel^borlng atoms must take place in order to achieve a stable con- 
figuration. Figure 5(b) shows the rearranged condition. Comparison 
of figure 5(b) with the positive dislocation in figure 4(b) shows 
that the two are essentially the same. The new positions of atoms 
5 and 6 after rearrangement were found by using the approximate 
method used by Taylor (reference 6) in his discussion of the motion 
of a dislocation under stress. The contributions to the potential 
field in row B by the neighboring rowo A and C etro takon to bo sinu- 
soidal with periods equal to the spacing of atoms 1 and 3 in row A 
and atoms 2 and 4 in row C. Figure 5(c) ^ows ■that addition of the 
two sinusoidal potential functions gives a resultant potential 
fxmetion having four minima. Atoms 5, 6, and 7 can bo oipoctod •to 
fall into the deepest of these minima if atom 5 had begun by oscil- 
lating to the l^t as in figure 5(a) . 

The displacement 5 of atoms 1 and 2 will be soms fraction f 
of 'the atomic separation dj^ in the slip direction. A value of f 
cf somewhat less than one -half makes it possible to form a stable 
ccxif iguratlon •with a minimum acti-vatlon energy. 

In the calculation of R„, it will be necessary to use Eyting's 
general theory of rate processes. A rl^t-hand source, under a 
shear stress Tg having -the direction shown in figure 5, will be 
considered. The "reactant" is the perfect lattice at the region of 
high stsr'ess ccxLcentration. The activated complex for a positive 
dlslocatlcvn, which for brevity ■will be called a positive activated 
complex, can be takon as ■the configuration shown in figure 5(a) and 
the product as the dislocation after atomic rearrangeanent has taken 
■place (fig. 5(b)). If ■the fundamen'bal assumption of rato-process 
theory that equilibrium exists between reactant and activated complex 
is made, the considerations involved in the general ■theory cf rate 
processes (reference 5) can be applied to the generation cf dislo- 
cations. From the theory cf rate processes, the rate of generation 
of positive dislocations can be expressed as 


1 
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The dlscuaslon la reference 5 (p. 189} shows that, If the dis- 
placement Involved in the transition from the activated ccmplex to 
the final state Is on the order of 10 centimeters, the factor l/t 
Is approxlmateljr equal to kT/h. In the fomatlon of a dislocation, 
this transition corresponds to passing frm the ccnf Iguratlon of 
figure 5(a) to that of figure 5(h), which cleeurly Involves a dis- 
placement of this order of magnitude. Hence, Is given approx- 

imately hy 

E+ = f P*- (8a) 

Similarly 

SJ = f P- (8b) 

The factors P"*" and P“ will now he calculated, lii order to 
simplify the calculation, the theory of elasticity Is assumed to he 
applicable on an atomic scale. The generation of a dislocation will 
he treated as a "local shear," which takes place h^r moans of thermal 
oscillations. The energy involved in this process will he calculated 
hy the theory of elasticity. Furthermore, Hooke's law will he assumed 
valid for large strains. The use of these aaBUS^tioas is necessazy 
hecause the forces hetwoen Individual atoms have not hoen evaluated 
In tezms of physical constants hut It should he remembered that the 
approximations may he serious enoxi^pi to lead to large errors In tlie 
results . 


The thezmal atomic oscillations In figure 5(a) will he reg£u:^lsd 
as resulting In a local shearing of atoms 1, 6, end 2 producing a 
strain of S/dg. Inasmuch as the elastic energy per unit volume for 
a shear stress T is t^/SG- and the volume of the region of shear 
is 2V, the elastic energy involved Is 


W = 


G 


( 9 ) 


The result will he applied to the generation cf a dislocation 
as follows: Let Tq he the local shear stress attained hy atoms 1, 

6, and 2 In figure 5(a) . If the displacement 8 is fd]^, appli- 
cation of Hooke's law to the dlsplEu^ement yields 

Gfdi 

To = -^=G3cf (10) 

^2 

The stress to he produced hy thermal motion Is 


•^t ' *^0 ' Ts 
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If it Is assumed tha-t the local shear stress to he attained h 7 
atans 1, 6, and 2 in figure 5(a) must be Tq, Independent of applied 
stress. The. stress depends csi the externally applied stress and 

Is In the same direction as Tq. Trcm equation (9) the energy to be 
attained thermally W.^. therefore Is 

= I (To - Ta)^ (11) 

Inaspiuch as Tg can be expected to be only a small fractl-m of Tq 
for rigid materials, 'the teim In Tg^ may be neglected relative to 
the other terms In the blnci&lal expemslon of (tq Tg)^. The 
approximation obtained Is therefore 

W.J. . VGr^f^ - EVbcTTg (12) 

vhere equation (lO) has been used. 

In order to evaluate Tg, vork-hardening must be oonsldered. 

A crystal in the steady -stato range of creep Is generally work- 
hardened because deformation has taken place. From the nature of a 
dislocation It is clear that a stress field exists about It. By 
making use of this stress field, Taylor (reference 6) has explained 
work-hardening oa the basis of a “lattice” of stuck dislocations. 
Taylor's explanation 1s a static one and Is applied only to deter- 
mining the form of the stress-strain curve. In terms of a dynamic 
approach, the existence of a stress field about a dislocation leads 
to the concept suggested by Eochendorfer (reference 11) that the 
lattice of stuck dislocations creates an Inner stress field whose 
direction at the point of generation is opposite to the externally 
applied shear stress Tg In the slip plane. The value of the inner 
stress field at the point of generation will be called the back 
stress T^. The dynamic concept of hardening therefore consists In 
the foimaticaa of a lattice of stuck dislocations, fdilch results In 
a back stress at the point of generation. This beick stress lowers 
the effective sheer stress at these points and as a result the rate 
of generation, hence the creep rate, decreases. These considerations 
lead to the result that 


Tg = a(Tg - T^) (13) 

I nasmu ch as t^ ie eJ.ways less than Tg, the directions of Tg and 
Tg are the same. 

In a large assembly of -systems, which eire in thermal equilibrium 
with each other, the fraction of them that will have energies greater 
than some arbitrary value e la given by e"W^. In addition to 
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the energy reguirement, the ocourrenoe cf the couf Iguratloa of fig- 
ure 5(a) InvolTee a llmltatlcoi oa the dlrectloa of atans 1^ Z, 3, 4, 
and 5 to a particular oryetallographlc direction vlthln a very small 
solid angle. If the prohahlllty of oscillation in this direction 
Is p, the total probability of the configuration of figure 5(a) ^ 
that Is, the probability that a rl£pit-hand source vlll form a posi- 
tive activated complex Is given by 

p+ p5e-Wt/KP . p5g-(7Gx2f2 - 2VxfTg)/KT (14a) 

Similarly the probability of forming a negative activated com- 
plex (the configuration of fig. 5(d)) at a positive so»aroe is 

P- . p5e-(VGx2f2 + 2VifTB)/i3? (14^,) 

Inasmuch as is opposite to Tq In this case. The applied 

stress Ta therefore makes the generation of a positive dislocation 
mare probable than the generation of a negative dislocation at a 
right-hand source. 

If a negative dislocation is generated in the region under dis- 
cussion, It would move to the left and out of the crystal without 
producing deformation if It were the only dislocation present. The 
last positive dislocation generated from this source, however, became 
stuck after moving through part cf the crystal lattice; that is, all 
forces upon It were balanced. The generation of a negative disloca- 
tion at the point being considered will upset this equilibrium and 
attract the previous positive dislocation. As the two dislocations 
approach each other, the stress that each exerts on the other 
increases (from the equation for the stress field about a disloca- 
tion, reference 6) and a union tedees place resulting in the annihi- 
lation of the two dislocations and restoration of a perfect lattice. 
The net rate of generation which Is the rate of generation 

of positive dislocations that will eventually produce plastic defor- 
mtlOT, is Hg - Rg or, from equations (8a), (8b), (13), (14a), and 

®g “ Q-('VGi2f2 + P«t)/m? ^2qVxf(Tg - T^,)/kT3 (15) , 

where 

P* « - k logg p5 (16) 

The equation for the net rate of generation of negative dislocations 
at a negative source is clearly Identical with equation (15) . Equa- 
tion (15) therefore gives the expression for the factor R of equa- 
tion (6b). ® 
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According to this derivation. It Is now clear that the genera- 
tion of a dislocation 2 atoms In length involves the slmulteineous 
oscillation of twice as many atoms with sufficiently great energies. 
Thus, the prohahility of such a generation taking place can be neg- 
lected as compared with the probability of generating a 1-atom 
dislocation. 

The hyperbolic sine term in equation (15) etrose from the lower- 
ing of the potential-energy barrier In the direction of the applied 
stress. The heat of actlvatlcn per molecule (the hel^t of 

the barrier for zero applied stress) and the entropy of activation 
per molecule ASg^ are (by comparison with equations (1) and (2)) 

= TCi^f^ (17a) 

ASg^ » - P* (17b) 

P* la always positive because p is less than. 1 (equation (16)), 
therefore is negative. 

When the general espresslon for Rg given by equation (15) is 
combined with the creep equation (6b), the ccanpleto equation for the 
steady-state rate of creep becomes 

U = !^i S (o - STi,)/^] (10) 

Li ll 

where the applied tensile stress a has been substituted for 
assuming that slip in the individual grains occurs in the plane of 
maximum shear stress. For hig^ values of applied tensile stress O 
the Iiyperbollo sins function, to a good approximation, can be replaced 
by an exponential and equation (18) becomes, after tearing logarithms 
of both sides 

log, u. . lo^ f) - ^ ( 0 - 2T,) (19) , 

Equation (18), or its approximation equation (19), is the form of 
the theoretical equation to be used for the steady-state creep rate 
of pure polyorystalllne materials as a function of stress and 
temperature . 


COMPARISON OF TBEOBETICAL AND EXPERIMENTAL CREEP EOLATIONS 

The evaluation of the back stress Tij as a function of temper- 
ature and applied stress is difficult and has not yet been acoanpllshod. 
A great deal may be learned, however, from a cotbparlson of the theoret- 
ical equation (19) with empirical equations. It has been empirically 
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detezmlned tliat a hyperlsollo siae function. rep3?esents the dependence 
of creep rate on applied tensile et3raas a for many polycrystalline 
materials (reference 12) In agreement vith equation (18). At hlg^ 
values of stress, the hyperhollc sine furtotlon reduces to ein expo- 
nential. In reference 13 creep data are presented, vhlcb show that 
a relation of the form 

loge (?) “ - I “ ■b + 70 (20) 

represents the experimental data in the range of hl^ stresses, 
where a and h are ccxistants for a given material and y Is a 
fimctlon of temperature hut not of stress. If the notation is so 
changed as to make this equation similar to equation (19), equa- 
tion (20) can he rewritten 

/^*^i kjr'N AjtBT ^ ea 

where 
A B ka 

B a kh + k lofi^ 

P = yKT 

Comparison of the forms of equatlcxis (19) and (21) show good agree- 
ment . The simultaneous validity of the two equations requires that 
T|j he a Uneeu: function of stress. In addition, the oondltlan 
T-jj = 0 for a a 0 leads to the following expression for T-jj 

a F(T)a (22) 

I 

Equation (21) shows that three factors completely determine the 
steady-state rate of creep; the constants A and B and the tem- 
perature function p. A comparison, of equation (21) with equa- 
tions (19) and (22) shows that these factors are expressed in terms 
of x^sical con^ants of the material as 

A - VGq f^x2 (23) 

B = P» - Vr^xSQ^a (24) 

P = qVxf(l - 2F(T)) (25) 

where the modulus of rigidity G has been taken in tho approxlmato 
tana. 

G = - oT) (26) 
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COMPUTATIONS AND PEBCISION 

Creep data were obtained from, literature for a number of poly- 
crystalUne "pure" metals that had been annealed (references 13 
to 20) . (Althou^ the references do not generally state the Impur- 
ity content, the metals will be considered pure eis compared with 
alloys.) The reproducibility of the creep data from the various 
sources is not known. From these data, curves of log creep rate 
against stress were plotted and strai^t lines were obtained. In 
tezms of equation (21), the coefficient 0/kT divided by 2.3 Is 
the slope of this straight line and the intercept Is equal to 

® VL h.J 2.3 kT • 

The values of ^ and E (defined as A BT) were then deter- 
mined from the slopos ».rid the Intercepts, respectively. Finally, 

A and B were separately determined from a linear plot of E 
against T whore d^^/L was taken as 10“^ for all materials. 

The oiTors Involved In the deturml natlcai of A, B, and 5 
depend In paart on the number of strosses for which data wore avail- 
able at each temporature and on the number of temperatures at which 
obsorvatlons were made. The ranges in temperatures and strosses 
used will also affect the error. Generally, the orrors for A, B, 
and 3 were estimated to be about tS x 10”13 o-rgs, t8 X 10"^6 ergs 
per 9K, and t30 percent, respectively. 

For a largo number of pure polycrystalline mstals. It v&b found 
that 3 varied with tomperatuz*e . The function 

3 = Cj^e®2^ 

approximates the temperature dopondence of 3> This function was 
determined over canporatlvely short ranges of tcanperature (200° to 
300° E) and therefore may not apply over a large temperature range. 

Values of the modulus of rigidity G wore obtained frem irefor- 
ences 21 and 22. These references shovod that variations in the 
value of G for a given material, ore about t5 i>ercont. Boom- 
temperature values wore used as approximations of the absolute-zero 
values of G. 


RESULTS AND DISCUSSION 


Factor A 

Equation (23) gives the theoretical expression 

VG^2j2 
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for the factor A of the emplr5.cal eqwtion (21) . In figure 6, 

A/V Is plotted as a function cf Gq The line drawn Is the hest 
strai^t line throu^ the origin and represents the data within the 
experimental error of A. From the preTlous theoretical discussion, 
the slope of this line is f^ and is equal to 0.25 or less. Because 
the activation energy predicted on the hasls of the validity of 
Hooke's law will be greater than the actual value. It can be expected 
that the theoretical prediction will be too hl^. The experimental 
slope from the sbral^t line of figure 6 is 0.14. 

The agreement between theoretical pr>jdlotlon and experimental 
results Is encouraging when the assumptions made and the variety of 
sources of data a:re considered. 


Factor B 

Vhen equations (19) and (22) wore compared with the empirical 
equation (21), B was found to be equal to 

pt - Tx2f2QQ0t. 

where a la determined frem the approximate relation 


G = Gq( 1 - oT) 

The order of magnitude of the term Vx^f^GQCt, which is 3 x 10"^® ergs 
per ?El, is within tlie estimated error of B, which is +8 x 10"^® ergs 
per this term can therefore be neglected In comparing the 

B values for different materials. The values •£ B for different 
meteJ-s are presented In the following table: 


Material 

B 

(ergs/°K) 

Beforence 

Aluminum 

4:0 X 10~^® 

13 

Copper 

45 ■ 

14, 15 

Iron 

39 

17 

Lead 

54 

20 

Nickel 

33 

16 

Platinum 

42 

13 

Silver 

36 

13 

Tin 

47 

19 


A correlation exists between tlie scatter of data presented In 
figure 6 and the scatter cf B values In this table. 03io straleht 
line in figure 6 falls between the oxperlmental points for aluminum 
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and iron. Tlie table sbowB that aluminum and iron hare B valuee 
equal to 40 x 10"^® and 39 x 10"^® ergs per respeetively . The 
points above the line In figure 6 are all related to metals that have 
B values leas than 40 X 10"^® ergs per 9 e and the points bej.cw the 
line correspond to B values greater than 40 x 10"^® ergs per °K. 

The correlation between the direction of scatter of data for A 
and B Is satisfactorily explained by considering the method used to 
obtain experimental values of A euid B. For each material, values 
of E (defined as A 4- BT) were plotted against T and the values 
of A and B were obtained from the intercept and slope, respec- 
tively, of the best straight line drawn throo^ the points. For a 
given material, the value -'of the intercept A viU increase as the 
slope B of the line drawn through the experimental points is 
decreased. In view of this explanation, the scatter of data for A 
about the theoretical line appears to result from experlmenbal error 
alone. It can also be concluded tliat tho correct values of B fall 
about their mean value much more closoly than those in the table, 

Vhen the term containing a is neglected, B becomes equal to 
P*, which from equation (17b) is the negative of the entropy of 
activation AS^. Ihe hl^ negative value of the entropy of activa- 
tion In the creep process has been pointed out In references 3 and 13 
but satisfactory q,uantltatlve explanations of the magnitude of this 
factor could not be given. 

On the basis of l^e present treaimiont, the restriction of dlrec- 
tlcHi of 5 atoDie participating in the generation of an activated com- 
plex accounts for the large negative entropy of activation. A good 
approximation of the value of the probability of oscillation p of 
on atom in llie slip direction can be obtained by substituting values 
of B for P’ in equation (16). A value of about 3 x 10"^ was 
obtained for p when P* was equal to 40 x 10"^® ergs per 


Factor ^ 

The theoretical expression for the factor 6 is given by 

P = qVXf(l - 2P(T)) 

The term F(T) in this equation is related to the back stress 
(equation (22)), which resulted from the inner stress field about a 
dislocation. This stress field is rolated to the modulus of rigid- 
ity Or (reference 6). Therefore, p at any given temperature was 
expected to bo a function of G, Figure 7 shows that at room tenper- 
aturo p varies as whore n is approximately 3. 
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In. order to detezmlne Aether tlie temperature dependence of G 
could alone account for the temperature dependence of 3, the varl- 
ation cf G vlth temperato^ vas InTestlgated. It can he seen In 
figure 8 that fcac Iron, 3 and G'® have about the aoms tenq>erature 
dependence. Figure 6 ahowa, however, that 3 la proportional to 
G'3 . Therefore, it can bo ocaicluded that F(t) contalna a temper- 
ature-dependent feujtor In addltlca to G. 


EVAUJATION OF KESILTS 

The atatement la made In reference 9 (p. 344) that "a full 
Interpretation cf creep result a In tems of dislocation theory would 
be welccme but has not yet been achieved.” In the present report 
two contributions are made to a fuller understanding of creep in 
toims of both rate-process theory and dislocation theory. The first 
contribution la the quantitative dependence cf creep on physical 
c exist ants of the material; the second is an understanding of the 
origin of the large negative entropy of act.lvatlon ASq^, which was 
a source of difficulty In previous treatments (references 3 and 13 ). 

Althou^ the previous theories have been successful in showing 
the dependence of creep on stress and temperature, none has been able 
to show the quantitative dependence of creep on physical constants cf 
materials. Fixm a prewtioal point cf view, equations that relate’ 
creep to stress and temperature have useful engineering application 
In predicting the creep behavior cf a given material. These equa- 
tions can yield little useful information, however, to the metallur- 
gist who must find criteria to use In synthesizing alloys for bettor 
creep resistance. A kncjwledge of the dependence cf creep rate on 
physical constants cf material would yield such criteria there- 
fore have a practical as well as an academic interest. 

(i^e empirical equation that gives the dependence of steady-state 
rate cf creep on stress o and absolute temperature T can be taken 
as 

Xogu. 106 0^^)- 1^*5^ ( 21 a) 

The theoretloal treatment in this report has shown that the param- 
eters A, B, and 3 depend on ihysioal constants of the material 
In the following manner: 


A = VGq f 2 x2 

B = P' - « P« - A50 

3 - qVSf(l - 2F(T)) ) 


(23) 

(24) 

(25) 
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Equation (21a) sIigwb "Uiat, in order to have a low rate of oreep, 

A B must he large and 0 must he small. Equation (24) 

Indicates that the effect of an? physical cccistant on A will he 
in the opposite direction to its effect on B. Ineismuch as d in 
equation (24) is c£ the order of 10“^ and P' heis heon ahovm to 
he approximately constant for all pure annealed materials, the phys- 
ical constant will have a greater effect on A than on B up to a 
temperature of 10^° E. The most Important physical constant in 
equation (23) Is the modulus of rigidity at ahsoluto zero Gq Inas- 
much as the other constants do not vary much In order of magnitude 
among different materials. Empirically 3 was also found to decrease 
with increasing Cr. Pure metals having hl^ moduli of rigidity will 
therefore he creep resistant. 

The relation between modulus of rigidity G and modulus of 
elasticity E is given hy 

r » B 

® 2(1 4- u) 

Inasmuch as values of (i usually run from about 0.3 to 0.4, the 
effect of variations in this factor for different materials will he 
so small that in general, pure metals will he creep reslsteuit if 
their modulus of elasticity Is hl^. 

The extension of the theory to alloys. In particular heat- 
resisting alloys, will require modification to include the effects 
of strains produced hy solid soluble atoms and precipitated peart 1- 
cles. In view of the need for criteria that will hasten the devel- 
opment of alleys having heat-reslstiiag properties, any indications 
that can he obtained from the present theory should ho considered. 

The theory suggests the use of materials of hl^ moduli of rigidity, 
and therefore of high moduli of elasticity, such as tungsten, molyb- 
denum, and cemontod txmgsten carbide as matrix znatorlals for heat- 
resisting alloys. An investigation conducted under the auspices of 
the National Defense Research Committee has shown that chrome-base 
alloys having high tungsten or molybdenum content generally have 
better heat -resist j.ng properties than currentlj' used alloys. The 
use of still hl^er peroontagos of those elomontB therefore merits 
serious consideration. 


Aircraft Engine Research Laboratory, 

National Advisory Coomittee for Aeronautics, 
Cleveland, Ohio, January 16, 1946. 
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Figure 1.- Typical creep curve showing three stages. 
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Figure 4. - Generation and movement of a dislocation. In 
sketches a, b, c, and d a positive dislocation moves to 
the right; In e, f, g, and h a negative one moves to the 
left; the resulting deformation is identical. (Taylor.) 
(Reproduced from fig. I of reference 9, p. 335.) 
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Figure 6«> Generation of a dislocation at a positive source. 
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Flgur* 8.- Tenperkture dep*nd«ne« of f and l/O for iron. 
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